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ABSTRACT 

Aims. The motivation for this study is to include the effect of plasma flow in Alfven wave (AW) damping via phase mixing and to 
explore the observational implications. 

Methods. Our magnetohydrodynamic (MHD) simulations and analytical calculations show that, when a background flow is present, 
mathematical expressions for the AW damping via phase mixing are modified by the following substitution: C' A (x ) —> C A (x) + Vq(x), 
where C A and V 0 are AW phase and the flow speeds, and the prime denotes a derivative in the direction across the background magnetic 
field. 

Results. In uniform magnetic fields and over-dense plasma structures, where C A is smaller than in the surrounding plasma, the flow, 
which is confined to the structure and going in the same direction as the AW, reduces the effect of phase-mixing, because on the edges 
of the structure C' A and V' Q have opposite signs. Thus, the wave damps by means of slower phase-mixing compared to the case without 
the flow. This is the result of the co-directional flow that reduces the wave front stretching in the transverse direction. Conversely, the 
counter-directional flow increases the wave front stretching in the transverse direction, therefore making the phase-mixing-induced 
heating more effective. Although the result is generic and is applicable to different laboratory or astrophysical plasma systems, we 
apply our findings to addressing the question why over-dense solar coronal open magnetic field structures (OMFS) are cooler than 
the background plasma. Observations show that the over-dense OMFS (e.g. solar coronal polar plumes) are cooler than surrounding 
plasma and that, in these structures, Doppler line-broadening is consistent with bulk plasma motions, such as AW. 

Conclusions. If over-dense solar coronal OMFS are heated by AW damping via phase-mixing, we show that, co-directional with AW, 
plasma flow in them reduces the phase-mixing induced-heating, thus providing an explanation of why they appear cooler than the 
background. 
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1. Introduction 

A large amount of work has been dedicated to understand the 
role of Alfven wave (AW) damping in providing heating for 
laboratory and astrophysical plasmas, be it the solar corona in 
general, or, more particularly, its open magnetic field structures 
(OMFS). Observed AW flux is sufficient to heat the corona ( As- 
chwanden l2005h . However, Spitzer resistivity alone is insuffi¬ 
cient to dissipate AWs efficiently dTsiklauri et al.ll2003l) . The 
phase-mixing of harmon ic AWs was propos e d as a way to al¬ 
leviate this problem by iHevvaerts & Priestl ( 1983 ). In phase¬ 
mixing, harmonic AW amplitude damps in time as B A w(x, 0 00 
Qxp(-rjC' A (x) 2 t 3 k 2 /6), where symbols have their usual meaning 
and C' A (x) denotes an Alfven speed derivative in the density 
inhomoge neity direction that runs a cross the background mag¬ 
netic field (IHevvaerts & Priestl 19831) . The phase-mixing of AWs 
that have a Gaussian profile (as opposed to harmonic) along 
the background magnetic field has_slowerj>ower-law damping, 
B A w 00 U 3/ 2 , as established by iHood et aD ( 2002 ). A mathe¬ 
matically more elegant derivation of the latter scaling law has 
been provided bv lTsiklauri et al.1 ( 2003 ). The exponentially di¬ 
verging magnetic fi eld lines provide even fas t er damping B AW = 
exp(-Ai exp(A 2 Q) dSimilon & Sudani 1 19891: iDe Moortel et al] 


l200(J : ISmith et al.ll2007l) . lMalara et al. (120001) considered small 

amplitude AW packets in Wentzel-Kramers-Brillouin (WKB) 
approximation in the Arnold-Beltrami-Childress (ABC) mag¬ 


netic field. The latter can (for certain values of physical parame¬ 
ters) have exponentially diverging magnetic fields, thus also pro¬ 
vidi ng a superfast (exponent of exponent) AW dissipation. Tsik- 
lauri (120141) studied the dissipation of AW in ABC fields using 
3D magnetohydrodynamic (MHD) simulations (without a WKB 
restriction) and found that perturbation energy grows in time. 
This was attributed to a new instability, whose growth rate ap¬ 
pears to be dependent on the valu e of the resistivity and th e spa¬ 
tial scale (wavelength) of the AW. Nakari akov et al.l (1 19981) stud¬ 
ied the nonlinear coupling of MHD waves in a cold, compress¬ 
ible plasma with a smoothly inhomogeneous low-speed steady 
flow that was directed along the magnetic field in the phase¬ 
mixing context. Their main focus, however, was on the wave¬ 
mode coupling rather than a possibility that the flow reduces the 
effect of phase-mixing, which we consider here. 

The OMFS in the solar corona - and possibly in coronae of 
other similar stars - come in different forms. We refer to the 
"openness" of the magnetic field in a sense that the structure 
must be able to sustain a background flow. These can be chromo¬ 
spheric upflows induced by magnetic reconnection in long coro¬ 
nal loops (long enough for the magnetic field to be treated as 
uniform, to simplify our model) or coronal polar plumes. Gen¬ 
erally, a distinction is drawn between coronal ho les (Golub & 
Pasachoffl2QQ9l). plumes ([Deforest et alJI 19971 : iDel Zanna et al. 
1997; Raouafi et al. 20071) and more recently dark jets (I Young 


2015) 
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The solar coronal holes (CH) are regions of low-density and 
low-temperature (compared to the background) plasma, which 
are believed to be a sourc e fo r fast_£y 800 km s -1 ) solar wind, 
(see Chapter 4.9 in lAschwanden ( 200 5). CH temperature is typ¬ 
ically 0.8 - 0.9 MK compared to surrounding quiet corona that 
has a temperature of 0.9 - 1.2 MK. The boundaries of CH can be 
clearly seen in soft X-ray images, because of the absence of hot 
1.2-1.5 MK plasma in them, when compared to the background. 

The solar coronal polar plumes (CPPs) are radial, thin elon¬ 
gated structures that are visible in white light eclipse pho¬ 
tographs as enhancements of density (3-6 times denser than the 
backgroun d), usually located inside coronal holes ( Del Zanna 
et al. 1 19971) . Because they are denser means that the plumes ap¬ 
pear brighter than the surrounding media. In extreme ultravio¬ 
let (EUV) spectroheliograms they appear as shorter spikes near 
the polar limb. Some models and observations suggest that the 
plume plasma remains much slower and cooler than inter-plume 
plasma up to 2.0 R Q . Values above these plasma parameters start 
to approach the inter-plume values, matching them at about 
3.0 R q . The flow speed and temperature increase of the plasma 
inside plumes is sometimes observed dRaouafi et al.ll2007 ). The 
latter work explains the flow speed and temperature increase by 
a possibility of interaction of the CPP’s core with the faster and 
hotter inter-plume material. Generally, it is debatable whether 
CPP is the source o f fas t solar wind (see related discussion and 
references in [Deforest et al.l ( 1997 )). Ultimately, this is related to 
the question of whether CPPs have small dipolar magnetic field 
patches a t their base or are unipolar. Extensive work bv Defor- 
est et al. d 19971) suggests that CPP have unipolar magnetic fields. 
Their Figure 9, however, shows that, despite the magnetic field 
being unipolar, it is still patchy, or, and this is crucial for our 
model, the Alfven speed varies across the magnetic field, giving 
rise to the phase - mixin g of AWs. Also, on the theoretical side, 
lOfman & Davila! ( 1997 ) have shown that torsional AWs generate 
solitary waves non-linearly and these may play a crucial role in 
fast solar wind acceleration, which is a separate issue. 

The solar coronal dark jets (DJs) are relatively new features 
(lYoungll2015l) . The coronal jets have been known for some time 
to be a feature of solar coronal hole observations obtained in X- 
ray or EUV wavelengths. I Young ( 1201 5) shows examples of DJs 
that are essentially invisible in EUV image sequences but have 
a clear signature in Dop plergram s derived from an EUV emis¬ 
sion line. Interestinglv JCirtain et al.1 ( 2007 ) provide evidence for 
Alfven waves in solar X-ray jets. 

Chapter 8.2.3 in 1 Aschwanden ( 2005 ) provides an overview 
of observations of Alfven waves in OMFS. Because AW do 
not perturb density, they can be detected using spectral obser¬ 
vations. The non-thermal broadening of EUV coronal lines typi¬ 
cally shows bulk plasma speeds of 30 km s -1 (iDovle et all 19981) . 
For a typical AW phase speed of 1000 km s -1 , this means that 
AWs in OMFS have amplitudes of 3% of the background. There¬ 
fore AWs in OMFS are weakly non-linear. Different types of 
waves are present not only in the corona. iJess et all ( 2012 ) use 
high spatial, spectral, and temporal resolution images, obtained 
using both ground- and space-based instrumentation, to investi¬ 
gate the coupling between wave phenomena observed at different 
heights in the solar atmosphere. 

The motivation for this study is to include the effect of 
plasma flow in AW damping via phase-mixing and explore its 
observational implications. Section 2 provides the model and 
analytical calculations. Section 3 describes the two-dimensional 
(2D) MHD simulations that corroborate the analytical results. 
Section 4 lists the main conclusions and outlines suggestions for 



Fig. 2. Difference between background flow speed at time t as a func¬ 
tion of ^-coordinate, V 0 (x,y = y ma x/2, t), and its initial value at t = 0, 
Vo(v ,y ~ Jmax/2, 0), i.e. AV y = Vo(v, y — t) — Vo(v, y — ymaxl^i 0) 

for different time instances. Dashed curve corresponds to t = 5, dotted 
to t - 10 and solid to t - 20. This numerical run is considered for the 
fastest background flow, with D = 3 (as in panel (c) from Figure [TJ. It 
is clear that by t - 20 the flow speed difference is very small « 0.005, 
i.e. the flow stays intact and does not disintegrate. 


future observational validation of the theory that is formulated in 
this work. 


2. The model and analytical calculations 

We describe AW dynamics using MHD equations for cold (the 
background pressure is equal to zero, po = 0 ), incompressible 
(density perturbation is equal to zero, p' = 0 ) plasma with non¬ 
zero resistivity and zero viscosity (77 ± 0 , v = 0 ): 


dV 

— + (V-V)V = 


(WxB)xB 

P 0 P 0 


( 1 ) 


— = V x (V x B) - 7/V x (V x B). (2) 

dt 

The background plasma quantities are denoted with subscript 
0, while perturbation (AW) are denoted with primes as follows: 
V 0 = (0, Vo, 0), B 0 = (0, B 0 , 0), po = Po W, and V' = (0,0, V' z ), 
B' = (0,0, B' z ). The z-coordinate is assumed to be an ignorable 
direction, i.e. d/dz = 0. 

Linearly polarised AW is described by the z-component of 

EqsQ]and[21 

W'z.„W'z ((V x B') x Bq) z 

-7- + V0-7- = -, ( 3 ) 

at dy jUqPo 


dB' 

—l=(Vx(V'x B 0 )) z + (V x (y 0 x B')) z + t]V 2 B' z . (4) 

In Eq.©, the vector identity V x (V x B) = V(V • B) - V 2 B has 
been used with the divergence of the magnetic field being zero 
(V • B = 0). Therefore, the system of equations that describes 
AW dynamics (and dissipation) is as follows 

dV' dV' B 0 dB' 

-*7 +Vo 1T = —IT’ (5) 

Ot dy P 0 P 0 dy 
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(a) 



X 


(c) 
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(b) 



X 


Fig. 1. Alfven (C A (x), dashed curve) and 
background plasma flow (Vo(*X dotted curve) 
speeds as a function of x-coordinate (across 
the magnetic field). Solid curve shows the sum 
of the two C A (x) + V 0 (x). The different panels 
show cases of (a ) D = 1, flat total speed pro¬ 
file across v-coordinate (no phase-mixing), (b) 
D - 2, forward flow exceeding AW speed, (c) 
D - 3 stronger forward flow, further exceeding 
the AW speed, and (d) D - -1 backward flow. 




Fig. 3. Shaded surface plots of B z (x,y) at different times for the case without the flow D - 0. Panel (a) is for t - 5, (b) for t - 10 and (c) for 
t = 20. Panel (d) is time evolution of AW amplitude, normalised to its initial value, for the same case. The solid line corresponds to the asymptotic 
solution for large times, Eg. (122k at the strongest density gradient point x = (907/3000) x (2/r) = 1.8996. A more general analytical form Eg. (1211 
is plotted with dashed curve for the same x value (we actually plot 5,(1.8996, y)/a 0 )- Crosses and open diamonds are MHD numerical simulation 
results in the strongest density gradient point v = (907/3000) x (2/r) = 1.8996 and away from the gradient v = (1/3000) x (2/r) = 0.0021 (the 
first grid cell in v-direction), respectively, by tracing the maximum value of the Gaussian AW pulse. Dash-triple-dotted line corresponds to the 
asymptotic solution for large times, Eq.JTTk which is independent of v. A more general analytical form Eq.(l30l) is plotted with dash-dotted curve. 
It is also independent of x because at the peak of the pulse the value of exponent is unity. 
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Fig. 4. The same as in Figure[3]but for the case of D - 1. 



dB' dV' dB' 


( 6 ) 


Next, we transform the equations into the frame co-moving with 
AW plus background flow speeds with the following coordinates 
(x,y, t) —> (x,g, r): x = x, g = y - CA(x)t - Vo(x)t and r = st 
(with s <$c 1). The derivatives in the new coordinate system are 
as follows: 


d_ 

dx 


d_ 

dx 


- (C A (x) + V' 0 (x))t—. 


(7) 


d _ d 
8}’~W 


(8) 


d d d 

7 T = ~(C A (x) + Vo(x))— + £—. (9) 

at 0 % or 

Applying the transform to the linearized first order system of 
equations © and ©, and algebraically cancelling the terms that 
contain Vo yields 


- qw I + 4 ) 1 '- 


Bq d_B[ 

PoPo d £ ’ 


( 10 ) 


d_ 


d \ , dV' z 

-Ca(x)— + s—\B z — Bq—^t + 


% 


id-CCM+W),!) B’, 


(ID 


In EqQj]in the V 2 = d 2 /dx 2 + d 2 /dy 2 we only kept d 2 /dx 2 be¬ 
cause of the phase-mixing d 2 1 dx 2 » d 2 1 dy 2 . The next step is to 
apply operator - C A (x)d/d £ + sd/dr to Eq.JTTT): 

2 7 d d\i 

- + £ - ■ 


d d 

~ CA(x) di* e si 


av: 


B'=B 0 \-C A (x)— +s— —f 


% 


+r] i -c A (x)i- + ] (JL 


d 


St ) \ dx 


(C A (x) + V' {) (x))t^ 


dr) d% 

\ 2 

B' 


We then apply operator d/dt; to Eq.(10) 


d 


~C A (x)— + £— j - 


d \ dVl Bo d 2 B' 


d% dr 


PoPo df 


and substitute the latter into Eq. (fl2li . 

“ c -* w | +e l) ^ 




d 2 B : 


d% <9r / \dx 


PoPo df 
2 


+17 -Ca(x)— + s— — - (C A (x) + V' {) (x))t— 


% 


K 


( 12 ) 


(13) 


(14) 


Eq.(fl4l) is an equation for B' z and can be solved analytically us¬ 
ing simplifying assumptions in the asymptotic limit of large time 
?/ta » 1. Here t a is the Alfven time t a = LjC A (x), with L 
and C A (x) = Bo/(poPo(x)) 0 ' 5 being a typical lengthscale of the 
system and Alfven speed, respectively. Ignoring s 2 1 order 
terms, whilst retaining only t 2 » 1 order terms in the term pro¬ 
portional to r\, yields 


d 2 B' 

- 2 C a (x)e ' 




(15) 
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Fig. 5. The same as in Figure [3]but for the case of D - 2. 



In the above equation, the term CA(x) 2 d 2 B' z /dZ 2 algebraically 
cancels out_ 

As in lTsiklauri et al.l ( 2003 ). we now introduce an auxiliary 
quantity that has a physical meaning of slow diffusion time for 
an AW: 


tKC' a (x) + Vq(x)) 2 t 3 = tj(CAA + VpW ) 2 ? 3 
6 s 3 6 


(16) 


and the derivative 


d_ _ 2s 3 a 

ds rj(C' A (x) + Vq(x)) 2 t 2 dr 


(17) 


Using the new notation and after integration by Eq.dTSb re¬ 
duces to the diffusion equation 


= d 2 B[ 

ds d% 2 


(18) 


Following similar approach as in iTsiklauri et al.1 ( 200 3). 
Eq.(fl8l) can be integrated as 


Eq.(l20l) generalises the well known Heyvaerts & Priest’s solution 
for the case of shear flow, which is modified by the following 
substitution C' A (x) —> C' A (x) + Vq(x). 

For a Gaussian pulse of the following mathematical form, 
B z (£',t = 0) = croe - ^ 2/2cr \ its substitution into Eq.([l9l) gives a 
solution 


B z = 


x exp 


__ 

^1 + r,(C' A (x) + y'(x)) 2 r 3 /3<r 2 

[y - (C A (x) + Vp(x))?] 2 
2(c 2 + tj(C' a (x) + Vq(x)) 2 ? 3 /3) ’ 


( 21 ) 


which generalises the solutions that have been obtained before 
( Hood et al . 2002 : Tsiklauri et al.ll2003l) . Here, ao = 1/5 ^flncr. 
In the asymptotic limit of large times, t , Eq.(l2T1) implies that the 
amplitude of AW Gaussian pulse damps as 


B z = 1 [2ni](C' A (x) + y'(x)) 2 /3]" 1/2 r 3/2 . 


( 22 ) 


3. Two-dimensional MHD simulations 


B. 


= -L_ f 

2^/ns J-o « 


r*+oo 


I exp 

J—oo 

4s 


B z (?,t = 0)d?. 


(19) 


Let us substitute a harmonic wave initial condition B z (g', t = 
0) = exp(/k£') into Eq.(fl9l). The integration provides a solution 


B 7 


~V(C' A (x)+v' U)) 2 t 3 k 2 /6 -ik(y-C A (x)t-V 0 (x)t) 


( 20 ) 


The 2D MHD numerical simulations of this work employ 
Lare2d lArber et~aD ( 200 1) - a Lagrangian remap code that 
solves non-linear MHD equations. Lare2d is second-order accu¬ 
rate in space and time. The code is available for download from 
http://ccpforge.cse.rl.ac.uk/gf/proj ect/lare2d/. 
Lare2d uses shock viscosity and gradient limiters to capture 
shock. However, the amplitudes considered in this work are 
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Fig. 6. The same as in Figure[3]but for the case of D = 3. The dotted line corresponds to the asymptotic solution for large times, Eq.d22l). at the 
strongest density gradient point v = (907/3000) x (2/r) = 1.8996, while solid line and dashed curve are kept the same as in Figure[3]for comparison. 


weakly non-linear. In all our numerical simulations we use a 2D 
box with 3000 x 24000 uniform grids in x and y direction, which 
have alength of 2n and 167T in each spatial direction, respec¬ 
tively. The distance, magnetic field, and density are normalised 
to their background values L, Z?o,Po- The velocity and time are 
normalised to the corresponding Alfven values C\ = 2?o/ y/PoPo 
and ta = L/C a at x = 0. Boundary conditions are periodic 
in both spatial directions. In all numerical runs, a normalised, 
uniform magnetic field, of strength unity, is in y-direction. The 
density has a profile in x-direction p(x) =1+9 exp(-(x - n) 4 ). 
Therefore, a normalised Alfven speed profile is 

C A {x) = 1/ Vpo(x) =1/^1 + 9exp(-(x - tt) 4 ). (23) 

Plasma beta and gravity are set to zero in all numerical runs. 
For numerical reasons, plasma beta is actually set to 10 -8 , but 
effectively it is zero. In all our simulations with AWs at t = 0 
we impose a Gaussian pulse which has two components, B z - 
0.01 exp( -(y-0 .5) 2 /(2x0.05 2 )), V z = -0.01 exp(-(y-0.5) 2 /(2x 
0.05 2 ))/ a/po(x), making it a linearly polarised AW packet with 
an amplitude of 0.01. The pulse starts at y = 0.5 and has a 
width of cr = 0.05. Only in Figure [7] the pulse starts at y = 8, 
so that the backflowing middle part stays within the simulation 
domain. Plasma viscosity is set to zero, while first and second 
shock visc osity coefficients are set at 0.01 and 0.05 (see Arber 
et al. (1200ll) for further details). The plasma resistivity is always 
set to rj = 5 x 10 -4 . The resistivity is quoted in units of p^LCa- 
Therefore, 1 /fj = S is the Lundquist number. Plasma flow runs 


along y-direction and its mathematical form is given by 

VoW = D-D/ = D- DC a (x), (24) 

with constant D = 0,1,2, 3 or D = -1 controlling the cases 
of (i) no flow (usual phase-mixing, ( Hood et al . 120021: Tsik- 
lauri et al. 120031) )"(ii) flat profile across v-coordinate (no phase 
mixing), (iii) forward flow, exceeding AW speed, (iv) stronger 
forward flow, further exceeding AW speed, and (v) backward 
flow, respectively. Each numerical run takes about five hours 
on 256 CPU computing cores on a COSMA-5 supercomputer 
http: //dirac. ac. uk/Resources. html The different exam¬ 
ples are illustrated in Figure [T] It follows from Figure [Ha) that 
for the chosen set of parameters, Alfven wave (Ca(x), dashed 
curve) in the over-dense region 1.5 < x < 4.5 lags with 
min(CA(x)) = 0.3162. For D - 1, the background plasma flow 
(Vo(y), dotted curve) speed has a maximum max(Vo(v)) =1- 
0.3162 = 0.6838 in the same region, such that the sum of the two 
Ca(x ) + VoM = 1 (solid curve) for all v. The same follows from 
the analytical expressions Eq.(l23l) and (l24b . i.e. for D - 1 the 
Ca(x ) terms cancel out. Thus, D = 1 case corresponds to a flat 
profile across v-coordinate (no phase-mixing). In other words 
the forward flow completely counteracts the wave-front stretch¬ 
ing because of the variation of Alfven speed with the transverse 
coordinate. Figure |TJb) is for D - 2, forward flow, whose speed 
exceeds AW speed reduction, such that max(Vo(v)) = 1.3675 
and max(CA(x) + Vo(x)) = 1.6838. The velocity sum differ¬ 
ence between over-dense and peripheral regions in Figure [Hb) 
is 1.6838 - 1 = 0.6838, the same as the Alfven speed decrease 
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Fig. 7. The same as in Figure[6]but for the case of D = -1. 

1 - 0.3162 = 0.6838. In other words, in Figure Qlb), the solid 
and dashed lines are perfectly symmetrical with respect to y = 1. 
Figure [He) is for D = 3, even stronger forward flow, such that 
max(Vo(x)) = 2.0513 and max(C A (x) + V 0 (x)) = 2.3675. Fig¬ 
ure Hid) is for D = -1 backward flow, such that min(VoM) = 
-0.6838 and min(C A (x) + VoW) = -0.3675. The velocity sum 
difference between over-dense and peripheral regions in Figure 
Etc) is 2.3675 - 1 = 1.3675, the same as that in Figure Eld) 
1 - (-0.3675) = 1.3675. In other words, in FigureEIc) and Fig¬ 
ure Eld), the solid curves have the same gradient strength, which 
is also clear visually. 

In Figure [2] we plot the difference between background flow 
speed at time t as a function of x-coordinate, Vo(x, y = y m ax/2 ,0 
and its initial value at t = 0, Vo(x,y = y max l 2,0), i.e. AV y = 
Vq(x, y = y max /2, 0 - V 0 (x, y = y max /2, 0) for different times. The 
dashed curve is for t = 5, while the dotted curve is for t = 10 
and solid for t = 20. In this numerical run the background flow 
is the fastest of all the performed numerical runs with D - 3 (as 
in panel (c) from Figure El) • We gather from Figure [2] that, by 
the end simulation time of t - 20, the flow speed difference is 
quite small, « 0.005. It is even smaller for earlier times and/or 
smaller values of D. Therefore in Figured we demonstrate that 
the background flow in the absence of the AW pulse does not 
break up. 

Figure [3] shows the numerical run results for the case of no 
background flow with D = 0. Thi s corresponds to a sim ilar 
set-up studied i nlHood et alJ (1 20021) : iTsiklauri et alJ (120031) or 
more recently in Tsiklauri ( 2014]), as the numerical code bench¬ 
marking exercise. Figure [3] shows shaded surface plots of the 



0.01 


AW, i.e. B z (x,y) at different times. Panel (a) is for t = 5, (b) 
for t — 10, and (c) for t — 20. Panel (d) shows the time evo¬ 
lution of AW amplitude, normalised to its initial value. The 
solid line corresponds to the asymptotic solution for large times, 
Eq.(l22l). A more general analytical form Eq.(l2ll) is plotted with 
a dashed curve. Crosses and open diamonds represent numer¬ 
ical simulation results in the strongest density gradient point 
x = (907/3000) x (2 n) - 1.8996 and away from the gradi¬ 
ent x = (1/3000) x (2 n) = 0.0021 (the first grid cell in re¬ 
direction), respectively. The Dash-triple-dotted line corresponds 
to the asymptotic solution for large times, Eq.(l3T1), which is in¬ 
dependent of x. A more general analytical form of Eq.(l30l) is 
plotted with a dash-dotted curve. It is also independent of x be¬ 
cause, at the peak of the pulse, the value of the exponent is unity. 
At t = 0 (snapshot not shown here) AW is initially flat as in, 
e.g. Figure ffia), but without the hump in the middle and instead 
it is located at y = 0.5, according to the initial conditions given 
above. The Alfven wavefront quickly damps (the shaded surface 
disappears from [3ja) to [3jc) in the density inhomogeneity re¬ 
gions x « 1.5 - 2.5 and x « 3.5 - 4.5, where the wave fronts 
distort strongly. The derivative of the Alfven speed, C' A (x ), in 
^-direction, which enters Eqs. (CHT) and (l22l) . with V 0 = 0 and 
Vq = 0 as D = 0, is responsible for the fast damping of the AW. 
Away from the density gradient regions, a much slower dissipa¬ 
tion takes place. The latter is hardly noticeable on the timescales 
concerned (0 < t < 20). Away from the density gradient regions, 
the analytical solutions Eqs. ([30 b and (fTfl) seem to match the 
corresponding numerical solution (open diamonds) well. 


Article number, page 7 ofITOl 






















A&A proofs: manuscript no. 27105 


In Figure[4]we present numerical run results in a similar man¬ 
ner to Figure© except that D = 1. Note that in paneled) we 
have adjusted the plot y-range to 0.4 - 1.4 to see the amplitude 
behaviour more clearly. At all times, we see that the AW front re¬ 
mains flat across x-coordinate, and the phase-mixing effect is ab¬ 
sent. Thus, for over-dense plasma structures, which have smaller 
Ca compared to the surrounding plasma, the plasma flow that is 
confined to this structure and running in the same direction as the 
AW, reduces the effect of phase-mixing because, on the edges of 
the structure, C\ and VJ have the opposite signs. In fact, we de¬ 
duce from Eqs.(l23l) and (l24l) that for D = 1, Ca(x) + Vo(x) = 
1 =const and therefore C A (x) + Vq(x) = 0. Thus, the rather slow 
wave damping is not due to phase-mixing but to the usual Spitzer 
resistivity (note the crosses nearly coincide with open diamonds 
in Figure 0Jd)). This small mismatch, of the order of 0.0008, 
could be due to the two following reasons: (i) AW has two com¬ 
ponents ( V Z ,B Z ) and, as stated above, the initial condition for V z 
is V z = -0.01 exp(-(y - 0.5) 2 /(2 x 0.05 2 ))/ ^p 0 (x). Thus, de- 
spite the fact that the wave front is flat, because of the presence 
of the flow, the density is still a function of transverse coordinate 
x; (ii) for the latter reason, the non-linearity will damp the wave 
front in a slightly different way. Although, our analytical prob¬ 
lem is linear, we validate it with a fully non-linear MHD code 
Lare2d. In Figured the maximal flow difference AV^ is approx¬ 
imately 0.005. Thus, the small mismatch of « 0.0008 could be 
attributed mostly to the non-linearity of the numerical code. We 
stress that the plasma resistivity is always set to f) = 5 x 10 -4 
everywhere. Consequently, the small difference is not due to the 
weak (logaritm of plasma parameter) density dependence of the 
Spitzer resistivity being affected by the background density pro¬ 
file. Thus, the wave damps only by Spitzer (uniform) resistivity 
not by phase-mixing. This is the result of the co-directional flow, 
which reduces the wave front stretching in the transverse direc¬ 
tion to zero, while the front remains flat across x at all times. 

Figure [5] shows numerical run results similar to those found 
in Figure [3 except for D = 2. The choice of value of D is such 
that wave front in Figure 0a)-[2c) now bends forward instead of 
backward compared to Figure Oa)-[3c). We notice that Figure 
(3d) seems identical to [3d). This is because the AW-amplitude 
damping, using the phase-mixing formula Eqs.(l22l). contains the 
square of the sum of the Alfven and flow speed derivatives. In 
the case of Figure [3 the wave front derivative is negative in the 
density inhomogeneity regions x « 1.5 - 2.5, while in the case of 
Figure [3 it is positive. However, since the derivative is squared, 
the net effect is the same - the case without a flow and with 
forward flow with D = 2, the AW damping is the same. 

In Figure [6] we present numerical run results similar to those 
found in Figure [3 except for D = 3. This is the strongest flow 
case we consider and it demonstrates that fast flow can induce 
wave front stretching to the extent that it exceeds the usual ef¬ 
fect of phase-mixing without a flow. The latter can be clearly 
seen by looking at the crosses and the dotted line in Figure [3d), 
which appear lower than the solid line, thus indicating a stronger 
damping, as prescribed by Eqs.(l2lT) and (123) . 

Figure 0 depicts the numerical run results as in Figured ex¬ 
cept for D = -1. This corresponds to the back-flow case, i.e. 
the AW propagates in the opposite direction to the background 
flow. Compared to the case without the flow, Figure [3a)-[3c), in 
Figure 0a)-0c) we see the AW stretching is stronger and there¬ 
fore the wave damping via the phase-mixing is faster. However, 
the damping is the same as in the case where D - 3. Therefore 
Figure 0d) appears identical to©d). 

To summarise, based on Eqs.(l23t and (l24k Ca(x) + VoM = 
D + (1 - D)Ca(x), and therefore, (C' A (x) + Vq(x)) 2 = (1 - 
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D) 2 C' a (x) 2 . The latter prescribes the AW damping via phase¬ 
mixing in Eqs. (l2TI) and (l22l) . Thus, (C' A (x) + F'(x)) 2 = 

C’Ax) 2 for D = 0, (Figure®, 

OrorD = 1, (Figure 0), 

C' A (x) 2 for D = 2, (Figured, 

4C^(x) 2 for D = 3, (Figure©, 

4C^(x) 2 forD = -1, (Figure©. 

Thus, the AW damping via phase-mixing is the same in the cases 
of D = 0 and D - 2, and similarly for cases D - 3 and 
D — — 1. For D — 1, the phase-mixing effect is zero. Note that 
phase-mixing always acts in addition to the usual (homogeneous 
plasma case) resistive damping, which diffuses the AW pulse in 
the other spatial direction (in y). This can be seen by the slight 
broadening of the pulse in y direction in Figures [30 

To quantify the AW Gaussian pulse damping in the homo¬ 
geneous plasma case, we now consider a one-dimensional (ID) 
analogue of our 2D model by suppressing the variation in x- 
direction. Therefore we consider a ID AW pulse moving along 
y in homogeneous plasma. In this case Equation (H is replaced 
by 




B^d 2 ^ 

PoPo d % 2 


+7 7 



a \ d 2 B[ 
+ £ dr) se ' 


(25) 


Here, (d/dx - ( C A (x ) + V' 0 (x))td / d^j operator has been re¬ 
placed by d 2 /dy 2 = d 2 /<9£ 2 because in the 2D phase-mixing case 
in Eq.(Q3]), we set d 2 /dx 2 » d 2 /dy 2 , while in ID case d/dx = 0. 
This means that, out of the V 2 operator, we retain d 2 /dy 2 . Fol¬ 
lowing a similar procedure, as described above, the equivalent 
form of Eq.([T5t is now 


-2Ca£- 


d 2 B : 


-riC A - 


d d 2 B' 7 


d^dr ‘ ^ d^ de 

We now introduce an auxiliary quantity, 

TjT Tjt 

Sl = 2e = 2’ 

Using the new notation and, after integration by Eq. 
duces to the diffusion equation 

9B[ d 2 B’ z 
ds\ de 
As above, Eq. 


(26) 


(27) 


re- 


(28) 


can be integrated using Eq.dT9l). For a Gaus¬ 
sian pulse of the following mathematical form, B z (e , t = 0) = 


ct 0 e 


b: 


,-r/2cr 2 


, its substitution into Eq.dT9b yields 


Qf 0 


Vl +2,? 1 /0- 2 


exp 


e 


2 (cr 2 + 2s i) 


(29) 


Substituting the definition of s\ provides the desired solutions 


B z = 


ct 0 


s]l + ptlcr 2 


exp 


tv - (Ca(x) + Vo(x))t] 2 


2 (cr 2 + pt) 


and its asymptotic limit of large times: 

B, = ; P^r 1 ' 2 r'« 


(30) 


(31) 


In the case of homogeneous plasma resistive damping of the 
Gaussian AW pulse, the amplitude damps as r 1/2 compared to 
r 3/2 for the 2D phase-mixing. In the uniform density regions, the 
analytical solutions Eqs. (l30t and (fTfl) seem to match the corre¬ 
sponding numerical solution (open diamonds) in all Figures [30 
well. 
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4. Conclusions 

This paper uses analytical calculations, corroborated by MHD 
simulations, to demonstrate that, when a flow is present, math¬ 
ematical expressions for the Alfven wave damping via phase¬ 
mixing are modified by the following substitution C' A (x) —> 
C' A (x) + Vq(x). In uniform magnetic field and over-dense plasma 
structures, in which Ca is smaller compared to the surround¬ 
ing plasma, the flow, which is confined to this structure, and 
in the same direction as the AW, reduces the effect of phase¬ 
mixing. This is because, on the edges of the structure, C A and 
Vq have opposite signs. As a result of this, the AW damping, 
via phase-mixing, is slower when compared to the case without 
the flow. For example, in the over-dense plasma structures with 
density inside ten times higher than outside, the co-directional 
with the wave flow with « 0 .7Ca in the middle of the over¬ 
density (see Figure [Ha) and[4]) can reduce the phase-mixing ef¬ 
fect to zero. This is the consequence of the co-directional flow 
reducing the wave front stretching in the transverse direction to 
zero. Conversely, the counter-directional flow increases the wave 
front stretching in the transverse direction, therefore making the 
phase-mixing effect more effective (see Figure [Hd) and [7]) com¬ 
pared to the case without the flow. The flows with « 1 .4Ca and 
« 2 Ca in the middle of the over-density, make the wave front 
go faster in the over-dense region compared to the surrounding 
plasma (see Figures[Hb)inc),[5]and[6l). In the case without the 
flow (FigureO, the wave front is slower in the over-dense region 
compared to the surrounding plasma. The dissipation of the AW 
via phase-mixing is: (i) the same for flows with « 1.4Ca as for 
those without the flow (although in both cases the wave front is 
bent forward and backward, respectively) and (ii) larger for case 
of ~ 2 Ca (which is probably unrealistic observationally, but is 
presented for completeness), as quantified by Eq.(l22l). We stress 
that the result is generic and is applicable to different labora¬ 
tory or astrophysical plasma systems where flows, density inho¬ 
mogeneity across the background magnetic field, and AW resis¬ 
tive dissipation are all present. Nonetheless, we apply our find¬ 
ings to address the question why over-dense solar coronal open 
magnetic field structures are cooler than the background plasma. 
Since observations show that the over-dense OMFS are cooler 
than the surrounding plasma, and that they are in regions where 
Doppler line-broadening is consistent with bulk plasma motions, 
e.g. the AW, we show that, if over-dense solar coronal OMFS 
are heated by AW damping via phase-mixing, the co-directional 
with the wave plasma flow in them reduces the phase-mixing- 
induced heating, thus providing an explanation for why they are 
cooler than the surrounding plasma. 

As mentioned in the introduction, and reiterated here, there 
is currently a disagreement as to whether CPPs are the source of 
fast solar wind (see related discussion and references in Defor- 
est et al. (119971) ). Some observations claim that the source of the 
fast solar wind is the inter-plume region. [Giordano et al.l (2 00d ) 
present a spectroscopic study of the ultraviolet coronal emission 
in a polar hole. They identify the inter-plume lanes and back¬ 
ground coronal hole regions as the channels in which the fast so¬ 
lar wind is preferentially accelerated, i.e. outside the plume the 
speed is higher than in the plume. We stress, however, that these 
observations only present three measurements of the flow speed: 
two on either side of the plume, where flow speed is found faster 
and one inside where flow speed is slower. This does not pre¬ 
clude a possibility that on the edges of the plume the flow speed 
falls to zero and thus Vq(x) > 0 inside the plume. Since, in the 
over-dense structures with a uniform magnetic field C A (x) < 0, 


the effect of the flow counteracting the phase-mixing is still a 
viable poss ibility. 

Harra et al. ( 2015 ) explored the changes in coronal non- 
thermal velocity (i.e. bulk flows or AWs) measurements at the 
poles from solar minimum to solar maximum, using Hinode 
EUV Imaging Spectrometer (EIS) data. They find that, although 
the intensity in the corona at the poles does tend to increase 
with the cycle, there are no significant changes in the V n0 n-thermai 
values. The locations of enhanced V n0 n-thermai values that they 
measure do not always have a counterpart in intensity, and they 
are sometimes located in weak emission regions. The next logi¬ 
cal step in corroborating our theory would be to check whether 
there is a correlation between temperature and non-thermal ve¬ 
locity in over- or under-dense OMFS. Care should be taken 
however in interpreting the observational results. We considered 
over-dense structures (e.g. coronal plumes) and found that, co- 
directional with the wave flows reduce phase-mixing and hence 
reduce heating/wave-dissipation, such that they appear cooler 
than their surroundings. Obviously, in under-dense structures, 
where AW phase speed inside the structure is higher, hence co- 
directional with the wave flows increase the phase mixing and 
thus increase wave-dissipation/temperature. This would result in 
a positive correlation between proton temperature and solar wind 
speed iTu & Marschl (1 9941) : (Horbury & Matteini, priv. comm., 
manuscript submitted for publication), i.e. faster streams are also 
hotter. Typically during periods of fast solar wind (V > 700 
km s -1 ), with fast streams (V = 150 - 200 km s _1 ), both the 
total magnetic field and the density are constant. Thus Alfven 
speed is constant across the streams. [Matteini et al.l (20151) show 
that in terms of the turbulence, the system appears to be in a lo¬ 
cal equilibrium, where there are no jumps in particle/field ener¬ 
gies. This behaviour suggests that the turbulence has evolved to 
a stage where the system is in equilibrium, making fast streams 
and background homogeneous. However, assuming that at the 
origin there were independent steams of plasma, with different 
physical properties, it is plausible that Alfven speed is different 
in different streams. It is likely that the condition \B\ = const 
dMatteini et al.l 120151) is the result of the relaxation of the tur¬ 
bulence and not the initial configuration, where jets could have 
had different Bo, which has subsequently been smoothed out. 
The next logical step would be to do a careful pressure-balance 
calculation, taking into account the temperature changes. 

Thus, based on our model, the ultimate factor for interpreting 
the observations is dependent on whether Alfven speed (which is 
a combination of both density and the magnetic field) is smaller 
or larger than in the surrounding plasma. In summary, when plot¬ 
ting a graph of coronal non-thermal velocity, V non -thermaU i.e. the 
background flow speed, versus temperature inside the structure, 
T , based on, e.g. EIS (Hinode)/ AIA (Solar Dynamics Observa¬ 
tory) observations in the solar corona or Helios observations in 
fast solar wind streams, the model predicts: 

a positive correlation of V non - t hermai with T in the case of 

structures in which Alfven speed is larger compared to the 

surrounding plasma; 

anti-correlation of V non -thermai with T in the case of structures 

in which Alfven speed is smaller compared to the surround¬ 
ing plasma (hence the title of this paper). 

These conclusions are based on natural assumptions that (i) AW 
phase-mixing has a major role to play in heating these structures 
and (ii) that the flow is forward (co-directional) with the AW (i.e. 
solar wind). There is also a caveat that in the above correlation, 
ynon-thermal means background flow rather than AW or turbulent 
motions, and that the disentangling of the two maybe difficult. 
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This is maybe easier on the solar disk rather than on the limb, 
because measuring Doppler shifts should enable us to differen¬ 
tiate between regular (up- or down-) flows from bulk turbulent 
motions, which only manifest themselves via the line broaden¬ 
ing. 
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